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Disordered systems have grown in importance in the past decades, with similar phenomena 
manifesting themselves in many different physical systems. Because of the difficulty of the topic, 
theoretical progress has mostly emerged from numerical studies or analytical approximations. Here, 
we provide an exact, analytical solution to the problem of uniform phase disorder in a system of 
identical scatterers arranged with varying separations along a line. Relying on a relationship with 
Legendre functions, we demonstrate a simple approach to computing statistics of the transmission 
probability (or the conductance, in the language of electronic transport), and its reciprocal (or the 
resistance). Our formalism also gives the probability distribution of the conductance, which reveals 
features missing from previous approaches to the problem. 



I. INTRODUCTION 

Disorder is ubiquitous in nature. The issue of universal 
properties of disordered systems has sparked much inter- 
est, both experimental and theoretical, beginning with 
Anderson locahzation in electronic systems, and by now 
including phenomena observable in a wide range of phys- 
ical systems. 

Transport in one-dimensional (ID) disordered systems, 
being the simplest problem, is naturally the first to be 
understood. Studies on electronic conductivity by mod- 
eling a ID wire as an array of potentials with random 
shapes and/or random spacings abound. For an excel- 
lent review of the state of affairs for ID systems in 1982, 
we refer the reader to the article by Erdos and Hcrndon.^ 
Overviews of subsequent developments can be found in 
Refs.HH andU 

Despite the comprehensive understanding of the ID 
problem, owing to the difhculty of the subject, exact 
analytical results are few and far between. Pertur- 
bative treatments, typically in the limit of weak scat- 
tering, are the main tools for analytical studies: the 
scaling theoryl^S^ the Born approximation for scattering,^ 
the DMPK equation,^ and many more. Exact results 
emerge mostly from numerical studies, or from the group- 
theoretic approaches of Refs. [U [51 and [31 which apply in 
general situations, but are nevertheless somewhat com- 
plicated because of the use of direct products of transfer 
matrices. 

In this article, we revisit the simplest problem of a 
ID single-scattering-channel system where the disorder in 
question is one of uniformly distributed phase disorder. A 
physical realization of a system manifesting such a disor- 
der comprises a stack of identical semi-transparent glass 
plates, with random spacings between adjacent plates. 
The task is to explore the transmission of a laser beam 
through the random stack. Because the typical separa- 
tion between plates is much larger than the wavelength 
of the light, a uniform distribution of separation will de- 
scribe the phase disorder well. Such a system was already 



considered by Stoked in 1862, who gave a ray-optical 
treatment. It was reconsidered by others, who applied 
wave optics, in particular by Berry and Klein^*^ in 1997 
(in fact, it was their work that inspired our current ef- 
forts). 



For this problem, an exact yet simple analytical treat- 
ment is actually possible, as we will present here. Disor- 
der average of all moments of the conductance and resis- 
tance (that is, the transmission probability and its recip- 
rocal) can be derived easily within a single framework. 
The same framework further permits direct derivation of 
the exact probability distributions for the conductance 
and resistance, which allows for computation of all statis- 
tics of the disordered system. This goes beyond past work 
that reconstructed the distribution for the conductance 
from its moments, ^ or derivations that relied on pertur- 
bative approaches. '^'^ The simplicity in our solution lies 
in a relation with Legendre functions, whose properties 
are well studied and are easily amenable to analytical 
manipulations. Our exact solution for this simple model 
can conceivably be used as the starting point for per- 
turbation towards other more realistic systems (see, for 
instance, the recent experimental proposal of Ref. [T2|l . 



The article will proceed as follows: After setting up the 
problem in Sectionjll] we describe (Section |lII[ ) how to av- 
erage over disorder using a recurrence relation, previously 
derived in Ref. [TBI This recurrence relation is applicable 
even for general disorder. Specializing to uniform phase 
disorder, in Sections|lV]and|Vj we make use of a close link 
to Legendre functions to write down closed-form expres- 
sions for the expected values of moments of conductance 
and resistance. Comparisons with existing results in the 
weak-scattering, long-chain limit are presented in Section 
VI The recurrence relation also gives the exact probabil- 
ity distributions for the conductance and the resistance, 
and these probability distributions are studied in detail 
in Section [Vll[ We close with a summary in Section [VlII| 
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II. PROBLEM SETUP: TRANSFER-MATRIX 
DESCRIPTION OF THE SCATTERING 
PROCESS 

Consider a ID chain of scatterers, such as a stack of 
semi-transparent glass plates, or a string of impurities. 
Both the scatterer and the scattered particle are assumed 
to be scalar particles, or at least only scalar degrees of 
freedom participate in the scattering process. The scat- 
terers are identical, but sit at locations with varying dis- 
tances between adjacent scatterers. The distances be- 
tween adjacent scatterers are the random variables in our 
system — the disorder. 

The scattering of the particle off the nth scatterer can 
be summarized by a transfer matrix T„ which propagates 
the wave function of the scattered particle past the nth 
scatterer (see Figure [T]) . T„ is obtained by solving the 
scattering problem for the nth scatterer, and accounts for 
multiple scattering off the same scatterer due to reflected 
waves from adjacent scatterers. For our problem, Tn can 
be decomposed into two parts. 



n-l 



Tn = D{kLn) T with D{ip) 



e"^ 








(1) 



where T describes the scattering due to the nth scatterer 
itself, and D{kLn) accounts for the phase acquired by the 
scattered particle in traveling distance L„ between the 
nth and (n-l-l)th scatterers. T takes the general form 



T = D{fi)T{{>)D{a) 

with rr,?^ - ( """"Ml"^) sinh(^ , 
with - y ^.^^j^^i^^ cosh(i^) 



2^) 



(2) 



Here, a and /3 denote overall phases, and T('d) contains 

[cosh(ii?)] and the re- 
12 = tanh(i?9). Since the 



the transmission amplitude t — 
flection amplitude r = \/l 
scatterers are identical, the same T describes every scat- 
terer. The disorder that stems from the variable separa- 
tion between adjacent scatterers is characterized by L„, 
for n = 1, . . . ,iV (with Ln = 0). 

The total transfer matrix, for a fixed configuration of 
N scatterers, is 



=TwTjv-i...ri 

'•■Pn-i 
2 



(3) 



where v?n = 2(a + /J + /cL„). T^^) can also be writ- 
ten in the form of Eq. ^ , with an effective A''-scatterer 
transmission amplitude i^^^ — [cosh(ii9'^^^)] ^, so that 



(4) 



with overall phases a^^^ and /S^^^. These phases and 
can be recursively expressed in terms of •& and the in- 
dividual spacings L„. This can be understood by adding 
one more scatterer to the end of the chain of N scatterers. 



n+\ 



N 



«71 + 1 



FIG. 1. The transfer matrix description of the scattering pro- 
cess. T is the transfer matrix for scattering due to a scatterer 
by itself (for example, a glass plate); r„ includes the effects 
of the spacing L„ between the nth and (n-l-l)th scatterers. 
^„{x) is the wave function of the scattered particle just be- 
fore it hits the nth scatterer, and can be written in terms of 
left- and right-moving waves as tpn{x) = ii„ 
where k is the momentum of the scattered particle. 



^e"="= -I- v„e 



The total transfer matrix is now r(^+i) = Tat+iT^^) = 
7^(/3(w+i))r(z9(^+i))i?(a(^+i)), with 1^(^+1) obeying 
the composition law 



coshi9(^+i) = cosh I? coshi?(^) 

-Hsinhz? sinhijW cos 



>N 



(5) 



Here, 4>n twice the phase sandwiched between the two 
matrices r(i?) and T{^'-^+^^), and is given in this case by 
(/>jv = 2{kLjsf -I- a -I- P^^^). We can regard 4>n, in place of 
Ln, as the parameter that represents the disorder in our 
system. Composition laws can also be written down for 
Q,(Af+i) g^jjfj ^(Af-i-i)^ they will not enter our analysis. 

As a convenient shorthand, we will use the notation 
C = coshf?, and S = sinhi? = y/C^ — 1. Similarly, one 
can define C^^) = cosht?(^) and S'(^) = sinhi?^^) so 
that the composition law appears succinctly as 



C(^+i) = CCW + cos(^^). 



(6) 



The total transmission probability (or dimensionless con- 
ductance, in the language of transport on ID wires) after 
N scatterers is 



TN 



c(^) + 1 



(7) 



The value of depends sensitively on the configuration 
of the scatterers, that is, on the values of the phases 0„. 



III. RECURRENCE RELATION: AVERAGING 
OVER DISORDER 

Rather than focusing on a particular configuration of 
the N scatterers, one is usually more interested in statis- 
tics of the entire ensemble of configurations. This re- 
quires a statement about the nature of the disorder, 
namely, the probability distribution d/i((/)„) for 0„. More 
generally, one can have correlated disorder, where the (f>n 
values are not statistically independent, but this is be- 
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yond the scope of our current discussion. The transmis- 
sion of N scatterers, averaged over the disorder, is then 



d^(02) / d^((/>i) 



(8) 

The dependences on (j)i,(j)2, ■ ■ ■ ,4>n-i are imphcit in 
C'(^). Writing out these dependences in full can be very 
complicated without being enlightening. 

To simplify the problem, we assume that the disorder 
is uniformly distributed, which permits the replacement 



2tt 



(9) 



in the disorder average. Here, the subscript (27r) de- 
notes integration over any 27r-interval of our choice. A 
uniform distribution is a good description whenever the 
separation between adjacent scatterers is itself uniformly 
distributed. It also describes the physical situation when- 
ever kL 3> 1 for some typical separation L between scat- 
terers. This often applies for monochromatic light scat- 
tered by a stack of semi-transparent glass plates with 
layers of air between them. In the case of a chain of ions 
trapped in a lattice, thermal motion of the ions within 
the trapping potential leads to random separation be- 
tween adjacent ions. The separation is usually concen- 
trated around some central value L, but for an electron 
incident with large momentum k, the phase fcL„ will ex- 
plore many cycles of 27r over even small deviations of L„ 
from L. A uniform distribution for (/>„ ~ fcL„ is then a 
fitting description. 



A. A recurrence relation 

To proceed with our analysis, we recall that the com- 
position law ^ allows us to compute (t^v) recursively 
with the aid of the recurrence relation derived in Ref. |T31 
We define a map Aic that transforms any function f{C') 
in accordance with 



{Mcf){C') 



(27r) 



2n 



f{CC' + SS' COS (f>). (10) 



Here, as before, C — cosh is such that t — [cosh(ii?)] ^ 
is the transmission amplitude for a single scatterer. C, 
the variable here, and S' = -\/C'^ — 1 are to be viewed 
as the hyperbolic cosine and sine of some Note that 
(Xc/)(C'-1)-/(C' = C). 



The map Mc possesses a symmetry that will prove 
useful later. For any two functions /(C) and g{C'), 



dC'{Mcf){C')g{C') 



dC'f{C'){Mc9){C'), 
(11) 



that is, we can consider Aic as acting on either / or g. 
To see this, observe that 



dC'iMcf)(.C')g{C') 

dC dC"f[C')]Cc{C',C")g{C") (12) 
1 Ji 



with the kernel 

K.c{C',C") = 



(27r; 



d(t) 

2tx 



8{GC' ^SS' cos cjj-C") 



■ [1 - - c''- 



C 



112 



(13) 



Here, 5( ) is the delta function, and the subscript + is 
such that [.t]^ = l/^/x when a: > 0, and is other- 
wise. Since K,c{C' ^C") is invariant under permutations 
of C, C" and C", it follows that the roles of / and g in 
Eq. ( 12 ) can be interchanged, as expressed by the sym- 



metry rule (11) 



From Ref. il3|, the average transmission probability is 

2 



(r^.) = (^^^-V)(C') 



for /(C) : 



(14) 



a fact that can be verified by writing out {tm) in full for 
N — 1,2,3,.... More generally, the disorder average of 
any function of C^^^ is given by 



/(cw)) = (x^-V)(c') 



C'=C 



(15) 



Another way of deriving Eq. ( 15 ) is to consider the 
probability density Wn{C') that C^^^ takes the value 
C , so that 



dC f{C')WN{C') 



(16) 



In view of the composition law ([6]), we have, for N 
m + with m, n positive integers. 



Wn{C')= r dC^Wrr^C^) r dC2Wr,{C2) f ^ 6 {CC^ + S^S^ COS cf> - C 



dCi / dC2Ty™(Ci)/Cc'(Ci,C2)W„(C2). 



(17) 
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m scatterers 



n scatterers 



FIG. 2. Illustration of concatenating two sub-chains of length 
m and n together to form a single chain of N = m + n scat- 
terers. 



The first line of Eq. (17) can be understood by splitting 
the N scatterers into two segments (see Figure [2]) , one 
comprising the left m scatterers, with an overall C^'"'* 
value equal to Ci, the other comprising the remaining n 
scatterers, with an overall C^"^ value equal to C2. The 
two segments are separated by a random phase and 
the composition law gives the overall C^^^ value for the 
TV scatterers. For Wi{C') 
iteratively W2{C'),W3iC') 



S{C' - C), Eq. ^ gives 
, summarized as 



Wn{C') = {McWn-i)(C') = {M^-'Wi){C'). (18) 



Upon inserting this formula into Eq. ( |16[ ) and recalling 
the symmetry Eq. (11) of Mc, we get Eq. (15) for the 
disorder average of /. 

As an example of the usefulness of Eq. ( [T5| , let us com- 
pute the average of logrAr. We set /(C") = log (prpi)- 
Performing the integration over cj) in {A4c'f){C'), we ob- 
serve that {Mcf){C') = f{C) + f{C'). Repeated appli- 
cations of A^c thus yields 



{\ogTN)=Nf{C) = N\ogTi, 



(19) 



where ti 



^q-j- is the transmission probability for 
a single scatterer. This expression comes as no surprise 
as log Tat is well known to be additive under disorder 
averaging, and in the current context, the result ( 
the main observation in the paper by Berry and K 



Note that formula ( 15 1 applies even for disorder that is 
not uniformly distributed, as long as we replace /(2^) 5^ 
in the definition of Mc by the more general / dfi{(j)). 



Equation (17) also holds for a general Wi{C") not neces- 



sarily equal to a delta function. 



B. Eigenfunctions of the recurrence map 

Computing the average of log r^r is easy because Mc 
acts in a simple way on the relevant f{C'). For more 
general functions, Mc can act in a complicated manner, 
and it becomes difficult to solve the recurrence relation 
directly to obtain a closed-form expression for the av- 
eraged quantity. The properties of the map Aic thus 
require thorough understanding before we can proceed 
further. Since we are to apply Aic repeatedly, eigen- 
functions of Mc — functions that remain invariant (apart 
from an overall factor) under the action of Mc — will be 
particularly useful. 



Consider the Legendre functions, A^(C"), which are 
functions that satisfy the Legendre differential equation 
(see, for example, Ref . fT4|) . 

^(1 - + '^(^ + 1) ^''(^') = 0- (20) 

For our current purposes, as the notation already sug- 
gests, C" is to be viewed as the hyperbolic cosine of some 
parameter, and we are interested in C between 1 and 
00. For a given there are two independent solutions to 
the Legendre equation, P,y{C') (Legendre functions of the 
first kind) and Qu{C') (Legendre functions of the second 
kind). Q^(C") blows up at C" = 1 while P^(C" = 1) = 1 
for all V. 

Suppose we begin with P^{C') and apply the recur- 
rence map Mc- We recall the addition theorem for 
Py{C') (see, for example, Ref. [T?|). 

P^{GC' + SS' cos (t>) (21) 

00 

= P,{C)P,{C') + 2 ^ P-'"(C<)P™(C>) cos(m0) 

m— 1 

where P™ are the associated Legendre functions, and 
C<(>) — min(max){C, C"}. Employing this formula in 
the integrand of {McPv){C) yields the eigenvalue equa- 
tion 



{McP.){C') = P,{C)P,{C'), 



(22) 



and the Legendre functions P^{C') are eigenfunctions of 
Mc- One can directly verify that {McPi/){C') satis- 
fies the Legendre equation with the same value of v, 
which permits writing {McP}^)iC) as a linear combi- 
nation of P^{C') and Q„{C'). Considering the value of 
{McPu){C') at C" = 1 leads to the conclusion (l22|. 



This eigenfunction property results in simple behavior 
of P^{C') under repeated applications of Mc, 

{M^,P,)iC') = P,{CrP,{C'). (23) 

The degree v can be any complex number. Of particu- 
lar importance to us are the cases when is a nonnegative 
integer, and when v takes the form v = — ^ + ix. When 
;/ = 0, 1, 2, . . ., we have the Legendre polynomials famil- 
iar from many areas of physics; the Legendre functions 
P_i+i2;(C") are known as the Mehler (or conical) func- 
tions, and have appeared in other physical problems, for 
example, the solution of Laplace's equation in toroidal 
coordinates. More relevant to our current subject, the 
Mehler functions were used in the exact solution of the 
DMPK equation (see Ref. |3] and references therein), as 
well a s in a different eigenvalue problem for the Anderson 
model.lSEll 



IV. MOMENTS OF THE RESISTANCE 

In the language of transport in electronic systems, the 
analogous quantity for transmission probability is the (di- 
mensionless) conductance. The reciprocal of the conduc- 
tance is the resistance, pN = ^/tn = h{C^^^ + 1) for 
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N scatterers. Our present concern is to compute the 
disorder-averaged resistance. We begin with f{C') = 
\{C' + 1) = \\Py(C') + Po(C')], which gives 



1 



(24) 



upon applying Eq. (23). For a long chain of scatterers, 
1, we have log(pjv) — A^logC, consistent with the 
usual expectation that the resistance grows exponentially 
as TV increases. 

For averages of moments of the resistance, (p^), for 
positive integer m, we start with /(C) = [^(C" + 1)]™. 
Noting that /(C) involves only positive integer powers 
of C", and every positive integer power of C" can be ex- 
panded in terms of the Legendre polynomials PeXC) (see, 
for example, Ref. fMlor lTS]) . we can similarly obtain (p^) 
without great effort by applying Eq. (p3). 

Here, f{C') - 



To illustrate, let us consider (p^) 



|(C' + 1)3 ^[P3(C") + 5P2(C') -f 9Pi(C') + bPoiC) 



From Eq. ( 23 ) , we then get 



20 



N 



N 



20 



Pi{C) 



N 



(25) 



an expression that would have been difficult to obtain had 
we tried to solve the recurrence relation directly. Note 
that the averages (p^) can be computed iteratively start- 
ing with m = 1 by employing recurrence formulas for 
Legendre polynomials that relate C'Pe{C') to Pg±i{C'). 

For large N, (p^) is dominated by the Legendre poly- 
nomial with the largest index £ = m, and hence. 



1^2 



N/ 



(ml) 
(2to)! 



P„(C)'^ for iV > L 



(26) 



Thus, log(p^) oc N for ^ 1 and any fixed m. Also, 
p^ has spread \/ (p^") — {p^n)"^, which, for large iV, is 
dominated by (p^")- 



A. The Mehler-Fock transformation 

The Mehler-Fock transformatiorP^ is an index trans- 
formation that makes use of the Mehler functions 
P_ 1 +i2:(C") as the basis for expansion. Formally, one 

defines the Mehler-Fock transform f{x) of f{C') as 

/•OO 

fix) = a:tanh(^x) / dC P_ /(C')- (27) 



This integral exists whenever f{C') is (weighted) square- 
integrablei^ such that 



ACWC -l\f{C')\ < OO. 



(28) 



The inverse transformation expresses the original func- 
tion in terms of its transform. 



/(C) 



dxP_ 



(29) 



The Mehler-Fock transform allows us to easily compute 
averages of physical quantities /(C*^^^), whenever / is 
square integrable in the sense of ([28]). We write, 

(/(CW)) ^ (X^-i/)(C') 



C'=C 



dx {M 



N-l 



dxP 



P-^+iJiC) fix) 
(C)^/(x), 



C'=C 

(30) 



where, exploiting the linearity of the map Aici we ap- 
plied A4c to the Mehler functions using (23). This gives 
a closed-form expression for (/(CW)), with two remain- 
ing items to evaluate — the Mehler-Fock transform of /, 
and the final x integral — but otherwise takes care of the 
complicated recursive composition law for concatenating 
N scatterers. 



B. Moments of the conductance 



V. MOMENTS OF THE CONDUCTANCE 

To compute the average transmission probability (or 
conductance), we note that the above method of expand- 
ing powers of C" in terms of Legendre polynomials no 
longer works, since the C" dependence in rjv occurs in 
the denominator. Whereas the Legendre polynomials are 
complete for — 1 < C" < 1, this completeness is of no help 
for C" > 1, which is the range of interest in the current 
context. 

Nevertheless, the idea of expanding in terms of Legen- 
dre functions still works. For studying the conductance 
and its moments, we expand, not in terms of Legendre 
polynomials, but in terms of Legendre functions of the 
Mehler type P_i^i^{C'), through the Mehler-Fock trans- 
formation. 



For moments of the transmission probability, or cquiv- 
alently, moments of the conductance, the Mehler-Fock 
transform f{x) can be worked out explicitly. For (t)J'), 
we set 

/(C')=(^)'" for^>i. (31) 

To compute its Mehler-Fock transform, we note that 
P_ 1 -i-ij;(C") can be written in terms of the hypergeome- 

tric function,i4 P^(C') = {^)'" F{-:y, ^) 
for = ia;. Inserting this into the Mehler-Fock 

transform integral for / gives 

/(x) = 2a; tanh(7ra;) (32) 
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after making the substitution y 
a standard onep^ 



■ This integral is 



f{x) — 22:tanh(7rx) Gm{x) 



(33) 



with Gm(x) 



r(m- \ +\x) 



r(m) 

sn2 , 



(m- 1)2 



G„i_i(x), 



where r(z) is the famiUar Gamma function. For m = 1, 
/(x) = 27ra;tanh(7ra;) sech(7rx), an expression that can be 
verified against known integrals involving Mehler func- 
tions (see, for example, Ref. [T^. 

With this, we have the expression for (t]J') for m > 1/2 
(including noninteger m values). 



' N I 



POO 

2 / dx xtanh{TTx) Gmix) P_ 
Jo '■ 



.+,JCy\ (34) 



The remaining integral over x can be estimated in 
limiting cases (for example, the weak-scattering limit 
discussed below), or evaluated numerically. We note 
that the Mehler functions are standard special functions 
for which efficient methods of numerical evaluation are 
known (for instance, via integral representations) or even 
built into mathematical software. 



VI. A LONG CHAIN OF WEAK SCATTERERS 

To connect with previous work on this subject, let us 
examine the limit of weak scattering, where we take C 
close to 1. At the same time, we assume a long chain of 
scatterers, so that iV ^ 1. 

We consider (r^J*), given by Eq. (34). For fixed rn and 



C, the integrand in Eq. (34) is significant only for small 
X values. This can be understood from the following 
integral representation for the Mehler function (see, for 
example, Ref. [T5]) . 



P 



(C) 



V2 



du 



cos{xu) 
\/C — cosh u ' 



(35) 



where C = cosh?? as usual. Since x enters only in the 
cosine in Eq. (35), for fixed C, P_i+ij.(C) stays bounded 
as a function of x. Furthermore, one can check that the 



factor a;tanh(7ra;) G„i{x) in the integrand of Eq. (34 1 van- 
ishes for large x due to the suppression from the factori- 
als in Gm{x) as x grows. These observations justify the 
approximation that, for C close enough to 1, or equiv- 
alently, for i9 close enough to zero, one can expand the 
cosine in Eq. (35) about xu — 0, and keep only the low- 
order terms. 



P_ 



[C] 



(C)e- 



P 



X 



dx' 



-P_ 



h+ix 



(C) 



(36) 



In the second line, we have defined 6 > such that 



and approximated 1 — bx'^ ~ e ''^ . 

With this, (t™) becomes much simpler: 

{T^^)=2[P_^{C)]''IM, 

f°° 2 

with Xmix) — / dx xta,nh{'Kx) Gm{x) e~^^^ 



(37) 



(38) 



For Nb 3> 1, the Gaussian suppression in Im{x) allows 
one to approximate the integral by Taylor-expanding the 
integrand about a; = 0. This gives 



-^m {x^ 



y/2 






r(m) 



(39) 



with equality attained when Nb oo. Putting the pieces 
together, we have 



' N I 



r(m) 



H 2 



1 / TT \3/2 

2\Nb) 



N 



(40) 



valid in the limit of a long chain of weak scatterers. This 
expression validates the conjecture of Ref. [121 namely. 



(t2) 



l/iN-2) 



P_i(C) asiV^cx), (41) 



since the quantity T(ri) = J^,^^ if (C + S cosip)~^/'^ in 
Ref. [13] is an integral representation of P_ i (C) \^ 

We can further approximate b and P_i{C) for <C 1, 



and b — 



1 



1?^ + 0(^^ 



16 



(42) 



From these, P_i(C) ~ e which upon inserting into 



Eq. (40), gives 



'N I 



r(m) 



2\Nb) ^ 



(43) 



This expression is identical to that found in Ref. 6 for 
the conduction of electrons in a ID wire, once we make 
the identification L/l — Nb, where L is the length of 
the wire, and / is the mean-free-path of the electrons. 



The factor of 1/4 in the exponential in Eq. (43 1 is the 
famihar ratio of log((rAr)) to (logTAr) = Aflogri(see, for 
example, Ref. [2]). From the expressions for P_i (C) and b 
in the weak-scattering limit, one can compute corrections 
to the standard answer of 1/4 as a function of 

Observe from Eq. (40) that, under the approximations 



above, (r^J*) cx {tn}, and that the proportionality factor 
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depends only on m, but not on the scattering strength. 
We can derive this statement more directly, using the 
following identityfSl 



r(m) 



amix) Gi{x) 



with ajn{x) = Y[ ( 1 



k = l 



{2x) 



{2k - 1) 



(44) 



so that we can write (without any approximation) 

2 



' N I 



r(m) 



(45) 



X dxxtaiih{Trx)a„,{x)Gi{x)P^i+i^{C)^y 

The expression within the curly braces is exactly (tn), if 
we can replace amix) by 1. Now, amix) — 1 whenever 
a; <C 1, and this is a good approximation whenever the 
integral in the curly braces gets its dominant contribu- 
tion from small x values. From our analysis above for 
the limit of a long chain of weak scatterers, the Gaus- 
sian suppression from the e~^^^ factor guarantees that, 
indeed, the integrand is important only for x ^ I. 
More generally, P_ i_|_ia,(C)^ is significant for x <C 1 

only whenever (sinhi?/??)^ ^ l (see Appendix |a|. This 
condition reduces to Nb ^ 1 in the limit of weak scat- 
tering. We thus conclude that 



(rN) 



r(m- ^ 
r(m) 



(46) 



whenever (sinh??/t9)"2' ^ 1, which holds when either 
N OT -d (or both) are large. This applies beyond the 
limit of a long chain of weak scatterers. For example, 
one can numerically verify the excellent accuracy of the 
above approximation when i!) — 5 and N — 5, for which 
(sinh-fJ/i^)^/^ ~ 850, where neither are the scatterers 
weak, nor is the chain a long one. This proportionality 
between moments of the conductance is a specific exam- 
ple of a similar statement applicable for more general 
types of disorder (but valid only when N oo) previ- 
ously pointed out in Ref. [31 



VII. THE PROBABILITY DISTRIBUTIONS 
FOR CONDUCTANCE AND RESISTANCE 

Armed with all the positive integer moments (t^), one 
expects to be able to reconstruct the probability distribu- 
tion for the conductance for given N and C. In Ref.I3J the 
probability distribution for the conductance was obtained 
by writing down an expansion of the Fourier transform 
of the distribution in terms of (t^J*). Approaching the 
probability distribution from a different angle, one can 



make use of our recurrence relation to propagate the ini- 
tial single-scatterer distribution to the distribution for N 
scatterers. Our approach uncovers features absent from 
a reconstruction via moments. 



We already have an expression (Eq. (18)) for M^jv(C") 
which we repeat here, 



WNiC) = iM^~^Wi)iC') 

= iM^-'W,,c')iCn 



C"=C 



(47) 



where, in the second equality, Wi_c'iC") = (5(C" — C"), 
and we have made use of Eq. One can verify that 

WNiC) > for 1 < C" < oo [m fact, WNiC) = for 
C > cosh(iVt?)], and that J^dCWNiC) = 1. From 
Eq. ( |47[ ), we immediately have the probability distribu- 
tions for the conductance r' = , and the resistance 



C' + l 



p' = l/r', related by Jacobian transformations. 



dCWNiC) 
dr'^{M^.-'W,)iC) 



C' = ^-l 



dp'2iM^-'Wi)iC) 



C'=2p'-1 



(48) 
< r' < 1, 

1 < p' < oo. 



Observe that Eq. BtI) bears resemblance with the ex- 
pression for (/(C(^))) in Eq. if we set /(C") = 
Wi^c'iC')- It thus seems plausible that our previous 
method of the Mehler-Fock transform might aid us in 
simplifying the recurrence formula for WNiC). The 
Mehler-Fock transform integral (over C" , with C held 
constant) gives fix) = a;tanh(7ra;)P_ i^i^(C"). Forget- 
ting for the moment that /(C") is not square-integrable 
in the sense of ( 28 1 , we employ the inverse-transform for- 
mula to yield 

/>C30 

WNiC) = J da;a:tanh(^x)P_i+i,(C)^P_i+i,(C'). 

(49) 

For computing statistics of square-integrable functions 
of C", we can safely regard Eq. (49) as a true identity 



for distributions, since we recover Eq. (30 1 when we re- 
place WNiC) in J^dC fiC)WNiC) with the right- 
hand side of Eq. ( [49| . For computing the average of 
non-square-integrable functions of C", like the moments 
of the resistance, the validity of Eq. (49) is less clear. If 
we take the long-chain (s — Nb — > oojf weak-scattering 
limit (i? ^ 1) of the right-hand side of Eq. (49), we obtain 



dx X tanh(7rx) P_ 1 iC)P_i (C) 



N 



~ -j= — / du (50) 

2v27r Jcosh-i(C') Vcoshu-C 

Here, we have approximated P^ij^^^iC)^ ~ e~''/^e~''^ 
as before, and also made use of an integral representation 
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for P_i+i^(C') of the fornP 



(C) = — coth(7rx) 



cosh 



sinfxw) 
du ^ — 
'^(C) vcoshu- 



(51) 

The Hmiting expression in Eq. ( |50[ ) is identical to the dis- 
tribution for C^^^ found in Ref.llllbv solving the DMPK 
equation, and also in Ref. |S] for the resistance in the An- 
derson model. This lends credibility to Eq. (49), even for 
moments of the resistance. Furthermore, the fact that 
/(C) = [(C -|- l)/2]'" is not square-integrable is perhaps 
not troublesome because Wn{C') — for C > cosh(iV??) 



(although not apparent from Eq. (49)), so that the C" in- 
tegral in (/(C(^))) cuts off at C = cosh(iVi?) rather than 
extending to infinity. 

A different source of concern lies with how the singu- 
larity of the delta function in Wi{C') propagates as N 
grows. For iV = 1, the singularity occurs at C — C; for 
TV = 2, it occurs at C = 1, since W2{C') = /Cc(C,C"); 
for N — 3, the singularity is at C" = C again (see Ap- 
pendix |b]); for A'' = 4, it goes back to C" = 1. In fact, for 
any even A^, we expect Wn{C' — 1) to be particularly 
large (or even infinite), because, as explained in detail 
in Appendix [Bj there exists an infinite family of phase 
configurations that attain C" = 1 whenever N is even. 
Since Wn{C' ^ 1) ^ Wn-i{C' = C), this large value 
of Wn{C' = 1) for A^ even is inherited by Wn{C' = C) 
for A^ odd. This again suggests that the large value of 
Wn{C') occurs at values that are different for A^ even and 
N odd, putting into suspect the existence of an asymp- 
totic probability distribution for large A. Fortunately, 
using our exact expressions for W^jv(C"), one can show 
that beyond A^ = 4, Wn{C') becomes finite everywhere, 
and the singularity present in small A values goes away; 
see Appendix [B] for a detailed analysis. This justifies the 
validity of the reconstruction of the probability distribu- 
tion of the conductance via its moments in Ref. 3, which 
relies on a Fourier transform that comes with an "equal 
almost everywhere" condition that automatically gets rid 
of singularities with zero measure. 



VIII. SUMMARY 

We have shown how one can analytically derive the 
statistics of a ID chain of identical scatterers of any 
length, with uniform phase disorder. Making use of the 
fact that Legendre functions are eigenfunctions of the 
recurrence relation governing the chain statistics, we ob- 
tained disorder-averaged moments of the resistance by 
expanding in terms of Legendre polynomials; disorder- 
averaged moments of the conductance came from expand- 
ing in terms of the Mehler functions via the Mehler-Fock 
transformation. The probability distributions of the con- 
ductance and the resistance can also be written in terms 
of the recurrence relation, and we pointed out singulari- 
ties absent in existing derivations. These extra features, 
despite the fact that they play little role in the physics in 



the end, remind us of the importance of exact and analyt- 
ical solutions, which may be the only way of identifying 
subtle properties and verifying the validity of common 
assumptions. 
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Appendix A: Mehler functions for large x 

We are interested in an approximate form for the 
Mehler function P_i^i^{C) applicable for large x. Fol- 
lowing Ref. 21, one can expand the Mehler function for 
large x in terms of Bessel functions. 



^ 71—0 



(Al) 



where the first three of the coefficients v4„('i9) are 



-i? coth ?9 - 1 



^ , (3i?cothi?+ 1)2 -8(2 + ^92) 
and A2{^)) = ^ ^ ^. (A2) 

One can check that the n — term gives, by far, the 
largest contribution, so 



P_ 



AC) 



Mxd) for a; > 1. (A3) 



sinhi? 

Let us consider [P_ i+ij.(C)]^, in light of this approx- 
imate expression. Since Jq{x'&) < 1 for all xi? values, we 
focus on the factor (z9/ sinh'!9)^/2_ Pqj- -ff <^ we have 

(i9/sinhi9)^/2 ~ n - i97l2j ~ upon recall- 

ing that 6 ~ ji?2 when ?9 < 1. This implies, for < 1, 
that [P_i^i^{C)]^ is small for large x whenever Nb ^ 1. 

More generally, [P_i^i^{C)]^ is small for large x when- 
ever 

/ sinh??^ 



V ^ 



»1, 



(A4) 



which can happen either when A^ ^ 1 or when is large. 



Appendix B: Singularities in Wn{C") 

Here, we examine in detail the argument outlined in 
Section VII regarding singularities in Wn{C'). We be- 



gin by studying how Wn{C") behaves for small A^ near 
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C = \ ioY N even, and C = C for N odd. Since 
Wy{C') = (5(C" - C), it is singular at C" = C. For 
W-i(C') = {McW^){C'), we have 



W^2(C') = /Cc(C, C) = - [(C - 1)(2C2 - C - 1)1 
vr 

(Bl) 

which is singular at C" = 1. Near the singularity, we have 

VF2(C')-^ ^^(C- 1)^1/2 forC'^1. (B2) 
Try 25* 

For Wz(C'\ we first write it in terms of an elliptic inte- 
gral. 



where a>«) ^ max(min) { [(CC + S-S") - \\\(C'^ + 
5-2) _ {CC - S'S")],45'35"}. 7^ is the complete elliptic 
integral,'^^ 

K{k) = / da 



■y/l — K^(sinQ;)2 
1 ,^2 



7r/2 p (^-h^ 



(1_k2)1/4 



(B4) 



which is finite everywhere for < k < 1 except at k = 1, 
corresponding to the point where C" = C. Expanding 
K{k) about K — \ (see, for example, Ref. [T^ . we can 
extract the behavior near the singularity, 



Nevertheless, for iV even, we expect W^iC' = 1) to be 
particularly large. This can be understood in the follow- 
ing way: C*-^-* = C" = 1 corresponds to perfect trans- 
mission through all N scatterers. For iV even, every con- 
figuration of phases 4>i^ . . . ^ipN-i such that 4>n/2 — t^, 
and (j)N/2-j + <f>N/2+j = 27r, for j = 1,2,..., N/2 - 1, 
gives C" = 1. There are thus infinitely many perfect- 
transmission configurations, since only the sum of each 
pair is constrained, but not the individual phases (apart 
from the central one). Furthermore, C is stationary, 
as the (/)„s vary, at C" = 1. This stationary property, 
together with the infinite family of configurations attain- 
ing C" = 1, will give a large, possibly infinite, value of 
M^AT at C" = 1. This large value of Wn{C' = 1) for N 
even is inherited by Wiq{C' — C) for N odd because 
Wn{C' = 1) = Wn-i{C' = C). This difference between 
N even and N odd in the location of the large, possibly 
singular value of Wm puts the existence of an asymptotic 
probability distribution as — ^ oo under suspicion. 

Fortunately, one can show that beyond = 4, Wn{C') 
becomes finite everywhere, and the singularity present for 
small N values goes away. The points of suspicion are at 
C" = 1 for TV even and C = C for N odd. We note that 
we can write 



W„,+„(C' = 1) 



dC"W^(C")Wn{C"), (B6) 



as can be shown using the symmetry property Eq. (11) 
of A^c- Then, we have 



WdC = 1) = / dC"W^{C"f 



(B7) 



Wi(C')^--—^\n\C'-C\ for|C'-C|^0. (B5) 

W4(C") for C" = 1 is also singular, with a behavior iden- 
tical to that of W^{C') in the vicinity of C = C . Beyond 
= 3, it becomes difficult to obtain a closed-form ex- 
pression for Wn{C'). 



Since W3(C") is singular at one point C" = C, the vici- 
nity of that singularity dominates the integral, but the 
logarithmic singularity behavior from (B5|, after squar- 
ing and integrating, becomes finite. This immediately 
also tells us that WsiC = C) = Wq{C' = 1) is finite. 
Thereafter, for larger N, we have Wn{C') finite every- 
where. 
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